Dynamical signatures of the vulcanization transition 
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Dynamical properties of vulcanized polymer networks are addressed via a Rouse-type model that 
incorporates the effect of permanent random crosslinks. The incoherent intermediate scattering 
function is computed in the sol and gel phases, and at the vulcanization transition between them. 
At any nonzero crosslink density within the sol phase Kohlrausch relaxation is found. The critical 
point is signalled by divergence of the longest time-scale, and at this point the scattering function 
decays algebraically, whereas within the gel phase it acquires a time-persistent part identified with 
the gel fraction. 



PACS numbers: 64.60.Ht, 82.70.Gg, 61.43. Fs 

Vulcanized matter has been suggested as a model sub- 
stance for the transition from a viscous fluid to an amor- 
phous solid state, i.e. the glass transition. Experimen- 
tally one has access to the static and dynamic critical 
behavior, and theoretically a reliable semi-microscopic 
model is available for actually computing properties of 
the system in the amorphous state as well as close to the 
glass transition. So far, the approach based on statistical 
mechanics has mainly focused on static properties. Here 
we extend the analysis to dynamic quantities and discuss 
the motion of a tagged monomer. We thereby show that 
vulcanized systems display many of the dynamic charac- 
teristics of glassy systems, in particular a Kohlrausch law 
(generalized to nonzero wavenumbers) in the fluid phase 
and a divergence of the longest time-scale as the glass 
transition is approached. 

The gelation transition has been discussed very suc- 
cessfully in terms of a percolation transition . The per- 
colation threshold is identified with the critical crosslink 
concentration, and the average cluster size with the 
weight average molecular weight in the sol phase. The 
percolation picture has been confirmed by a statistical- 
mechanical approach [g-Q]. In addition, the latter ap- 
proach allows one to discuss structural and elastic prop- 
erties of the gel phase, which can not be obtained from 
a percolation model unless additional assumptions are 
made. The shortcomings of percolation theory become 
even more apparent if one is interested in dynamic quan- 
tities. There is obviously no dynamic percolation model 
that could account for the observed anomalous relaxation 
while approaching the gelation transition from the sol 
phase. 

Here we start from a (semi-)microscopic equation of 
motion, the Rouse model (5), and compute the incoher- 
ent intermediate scattering function as the vulcanization 
transition is approached from the liquid side. We only 
consider permanently formed crosslinks and model them 
by a harmonic potential. We furthermore assume that 
each polymer is equally likely to be crosslinked to each 
other polymer. We do not distinguish between gelation 
and vulcanization and use sol phase and fluid phase as 



synonyms, because, as we shall see, the long time dynam- 
ics is not affected by the molecular size and structure of 
the building blocks of the random network. Neglecting 
excluded volume interactions, we analyse the dynamics 
by making use of the theory of random graphs, created 
by Erdos and Renyi j(| . (For an informal account of ran- 
dom graph theory see Ref. ||.) Our main results are the 
following: 

• We observe three distinct phases with qualitatively 
different long-time dynamics. 

• If the number of crosslinks M is submacroscopic 
(i.e. o(N), where N is the number of macromolecular 
chains in the uncrosslinked system) then the long-time 
behavior of the autocorrelation is diffusive: S t (q) oc 
exp{-D q 2 i}. 

• If the crosslink concentration c := M/N is non-zero 
but smaller than the critical value c cr it for vulcanization 
then the time-decay of the autocorrelation is character- 
ized by a spectrum of relaxation times, reflecting the dif- 
fusive motion of a distribution of cluster sizes. The auto- 
correlation decays for large times as a stretched exponen- 
tial: St(q) oc exp{— (i/rq) 1 / 2 } with, however, a scaling of 
time with wavenumber that is specific to a diffusion pro- 
cess: Tg 1 oc D q 2 . This relaxation is the true asymptotic 
behavior, and should not be confused with a formally 
similar time decay in the Rouse model for uncrosslinked 
systems, which occurs for time-scales smaller than the 
Rouse time j5j. 

• As the vulcanization transition is approached, i.e. for 
c T c crit, the longest time-scale diverges like the inverse 
squared distance from the critical point, whereas the gen- 
eralized diffusion constant goes to zero linearly and, more 
generally, is proportional to the weight average molecular 
weight. 

• Precisely at c cr it the distribution of relaxation times 
falls off algebraically for large relaxation times, implying 
a (Do? 2 ^) 1 ^ 2 decay of the autocorrelation. 

• In the gel phase, i.e. for c > c cr ; t , the autocorrelation 
acquires a time-persistent part, which is identical to the 
gel fraction as obtained from static calculations. 

• The mean square displacement grows linearly with 
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time for all crosslink concentrations, with only a weak 
singularity at the vulcanization transition. Its long-time 
behavior is dominated by the "mobile" monomers on fi- 
nite clusters, even in the gel phase. 

We consider a system of N linear, identical, mono- 
disperse chains of arc-length L with force-free ends. 
Monomer s on chain i is characterized by its time- 
dependent position vector R t (i,s) (i — l,...,N and 
< s < L) in ei-dimensional space. The simplest dy- 
namics is purely relaxational: 

1 5H 

d t Rt(i,s) = -- +q t (i,s). (1) 

C 6Rt(i, s) 

The monomers' relaxation to the stationary state 
5H/SR = with rate £ is perturbed by thermal fluc- 
tuations, which are modeled as Gaussian white noise 
r\ with zero mean and variance s) r\^, (i', s')) = 

2 £-i S(t-t') 8 i<v S(s - s'), where v, fi = 1, . . . , d. We 
measure energies in units such that fceT is unity. If the 
Hamiltonian were to consist only of the Wiener term 




then the dynamical equation would reduce to the Rouse 
model. Here I denotes the persistence length of the in- 
dividual polymers. In the following we wish to treat a 
crosslinked melt or solution. We only consider the case of 
permanently formed crosslinks that constrain randomly 
chosen pairs of monomers, and we study the relaxation 
of monomers in the presence of a fixed crosslink configu- 
ration. A particular realization of M crosslinks is char- 
acterized by a set of indices {i e , s e ; i' e , Se}e£i specifying 
which monomer is connected to which other monomer. 
We adopt here the simplest distribution of crosslinks, 
i.e. the crosslinks are independent of each other and for 
each crosslink the pair (i e ',i' e ) of polymer indices is re- 
alized with equal probability. We leave unspecified the 
distribution of monomer labels (s e ;s' e ) because its form 
will not affect our results. Thus, for instance, the case of 
end-linking with unconstrained functionality is included. 
The crosslinks are modeled by a harmonic potential 

, M 

U ~2^ EW.'.l-CO) 2 , (3) 

e— 1 

where a > controls the distance between monomers 
within a crosslink. The Hamiltonian that we use is thus 
H := Hw + U . Hard ^-constraints for the crosslinks can 
be recovered from the harmonic potential in the limit 

The quantity of interest is the intermediate incoherent 
scattering function 

(exp{iq- (Rt +to (i,s)-R to (i,s))}). (4) 



The above average refers to the noise r\ and is taken with 
fixed given crosslink configuration. 

As the Hamiltonian is quadratic the Langevin equa- 
tion (]l|) is linear and R is a Gaussian Markov process, 
whose distribution is in the limit t Q — > oo characterized 
by 

1 ((R t+to (i,s)- R to (i, s)) (R t+to (if, s') - R to s'))) 

= - /*dre- Tr ( l , S ;i', S ')=:G t ( l , S ;i', S ')- (5) 
s Jo 

Here e _rr (i, s; i', s') is the matrix representation of the 
exponential of the linear operator T which stems from 
the deterministic part of Eq. (Q) and is defined by its 
action, viz., 

^° e=l 

X (f(ie,S e )- f(i' e ,s' e )), 

on functions / over the configuration space obeying Neu- 
mann boundary-conditions with respect to s. 

We are only concerned with the long-time behavior of 
the scattering function. Thus we need the asymptotics of 
G t as t — > oo, which (because T is non- negative definite 
by inspection) is governed by the eigenspace of T corre- 
sponding to zero eigenvalues. Moreover, Tf = implies 
(/, Tf) — 0, and this in turn implies (by virtue of Eq. (||) 
and the boundary conditions) that both df(i,s)/ds = 
for all i — 1,...,N and f(i e ,s e ) = f(i' e , s' e ) for all 
e = 1, . . . , M. Let us group the polymers into K clusters 
{Jsfk}fc—i, which are defined as maximal path- wise con- 
nected components. Then, according to the above rea- 
soning, the eigenspace in question consists of functions 
that are constant when restricted to any one cluster. We 
therefore have, as t — > oo, 

- 1 

G t (i, *;<',*') = D tJ2jf^(i)^ (i') + 0(l), (7) 
fe=i k 

with the bare diffusion constant Dq = \/(C,L). Here Nj. 
denotes the number of polymers in the /c-th cluster A4, 
whence J2k=i^ k = and S_\f k (i) = 1 if i € Nk and 
•Wf, (*) = otherwise. 

Due to the Gaussian nature of the process R we can 
express the scattering function in terms of G, and find 
for the long-time behavior of the scattering function 

^(g) = ^E^exp|-q 2 ^+O fc (l))|. (8) 

We have equipped the long-time corrections Ofc(l) with 
the cluster index k to indicate that we are not allowed to 
handle them as an overall pre-factor. 
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The above result (g) is not at all surprising. It merely 
confirms the notion that the slowest modes correspond to 
the center-of-mass diffusion of individual polymer clus- 
ters. Here Do/Nk may be interpreted as the diffusion 
constant of the fc-th cluster, and Nk/N as the probability 
for a given site to belong to the fc-th cluster. We remark 
that in accordance with this notion the strength of the 
crosslink potential does not enter in Eq. (^). Thus the 
result covers the limiting case, a — > 0, of hard crosslinks, 
too. 

The statistics of the polymer clusters' sizes is deter- 
mined by the crosslink distribution. Two different poly- 
mers are considered to be linked if there is at least one 
crosslink joining any two monomers of the respective 
polymers. In our model we specify the total number 
of crosslinks M so that the total number of polymer- 
polymer links varies in different crosslink realizations. 
These fluctuations are suppressed [|l0| in the macroscopic 
limit as there are 0(N 2 ) possibilities for M = O(N) 
crosslinks to form M different polymer-polymer links and 
only O(N) other possibilities. For any given number of 
polymer-polymer links, the statistics of the clusters' sizes 
can be taken over from Ref. Q by identifying polymer 
clusters with components of random graphs. 

As shown in Ref. [||, the polymer clusters exhibit the 
analogue of a percolation transition. Below the crosslink 
concentration c cr it = 1/2 there will be no clusters of 
macroscopic size and almost all polymers are members 
of tree clusters |llj , which are defined by the absence of 
loops. For c > \ there will be O(N) polymers in exactly 
one macroscopic cluster, which is not a tree, and almost 
all other polymers still belong to tree clusters of sub- 
macroscopic size Jl2[ . Let T n denote the average number 
of tree-clusters consisting of n polymers. Then, whatever 
the crosslink concentration c > may be, one has [fl3| 

lim L Tn = r n = ^ r !— (2ce- 2c )\ (9) 
iv->oo N 2c n\ v ' 

In the macroscopic limit, the average fraction of poly- 
mers in the giant cluster is given by JTif 

oo 

<Soc = 1-^™t„ (10) 

n=l 

which, for c < |, amounts to = p5[ . 

To work out the long-time behavior of the autocorre- 
lation, given the statistics of polymer clusters' sizes, we 
shall first make the simplifying assumption that we can 
ignore the possible dependence of the long-time correc- 
tion Ofc(l) on the index fc of the polymer cluster. This 
assumption can indeed be justified |Tc| ] for finite cluster 
sizes. From Eqs. (||) and (|^) we get for the scattering 
function, averaged over the distribution of crosslinks 

oo 

S t {q) = Soc + ^2 nT nexp{-D q 2 t/n} . (11) 

n=l 



This result holds for all c > with Soc, = for < c <i . 
In fact, one would expect S t (q) to be self-averaging p0[ . 
We now discuss four distinct regimes: 

• For a submacroscopic number of crosslinks there is 
only a submacroscopic number of clusters, which do not 
contribute in the thermodynamic limit. Hence the scat- 
tering function decays, to leading order, according to the 
Rouse model: S t (q) = exp{— D q 2 t}. 

• In the sol phase with a non-zero concentration of 
crosslinks, i.e. for < c < ^, the distribution of cluster 
sizes falls off exponentially for large n: 

nr n ~ - 1 n~ 3/2 exp{-n/i(e)}, (12) 

with h(c) := 2c— 1 — ln(2c) > 0. This exponential distri- 
bution of cluster sizes — and hence diffusive time-scales — 
gives rise to a stretched exponential decay of the scatter- 
ing function in the sol phase: 

St(q) ~ exp{-2(M c )A)9 2 i) 1/2 } (13) 

for D§q 2 t — > oo. The time-scale of the stretched expo- 
nential diverges like (c— |) as the vulcanization tran- 
sition is approached, whereas the generalized diffusion 
constant, 

r°° n^ 1 

goes to zero linearly as c f ^. 

• At the critical point for vulcanization, i.e. at c = 
\ , the distribution of cluster sizes falls off algebraically, 
nr„ oc n~ 3 / 2 . Consequently the autocorrelation shows 
algebraic decay in time: S t (q) oc (D^t) -1 / 2 . 

• In the gel phase, i.e. for c > ^, the scattering function 
acquires a time-persistent part Soo , which is given fl5|| by 
the largest root of 1 — = exp ( — 2c S^) . As expected, 
the time-persistent part of the autocorrelation is just the 
fraction of polymers in the macroscopic cluster, i.e. the 
gel fraction. The distribution ([)[) of finite cluster sizes 
is valid on both sides of the transition ]lq] . Hence we 
expect to see the stretched exponential behaviour dl3| ) 
for the decaying part of the scattering function in the gel 
phase as well. 

Lastly, we mention that the time-delayed square 
monomer displacement, 

1 N r L 

Cf-=jfYl / dsG t (i,s;i,s)~2Dt, (15) 
i=i ~' 

increases asymptotically linearly with time for all 
crosslink concentrations c. In terms of polymer clus- 
ters the proportionality constant D is given by D = 
DqK/N — DoJ2^Li T n- This changes smoothly with 
crosslink concentration: in the sol phase it is given by 
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D = Dq(1 — c); at c = 1/2 it is twice continuously differ- 
entiable; deep in the gel phase we find D ~ Do exp(— 2c). 
Hence the squared monomer displacement is unsuitable 
for detecting the sol-gel transition because its long-time 
behavior is always dominated by the finite polymer clus- 
ters. 

We conclude with a few remarks concerning the exper- 
imental situation. Relaxation following a stretched ex- 
ponential exp{— (i/ T ) } seems to be a universal feature 
of glassy systems. It has been observed in fluids under- 
going either thermal or chemical vitrification, as well as 
in magnetic systems undergoing a spin-glass transition. 
It can always be interpreted in terms of a distribution 
of time-scales, although the nature of the underlying el- 
ementary process may be quite different. In the present 
model the elementary dynamical process is diffusion, and 
the stretched exponential results from a superposition of 
many diffusive modes, reflecting the distribution of clus- 
ter sizes [|l7| . This gives rise to a highly characteristic q- 
dependence of the time-scale of the stretched exponential, 
namely t^ 1 oc D$q 2 . Such a time-wavenumber scaling 
has indeed been seen in neutron spin-echo experiments 
on inorganic and organic glasses . 

As the gelation transition is approached one observes 
experimentally a divergence of the longest time-scale, and 
at the critical point density fluctuations are found to de- 
cay algebraically. We compare our results to the data of 
Martin et al. JlS( ] . These authors use quasi-elastic light 
scattering to study silica sol-gels in the pre- and post-gel 
regimes. They observe a stretched exponential in the sol 
phase with an exponent (3 w 0.66 and, furthermore, relate 
the time-scale r of the stretched exponential to a diffu- 
sive time scale. The latter is found to diverge as the gel 
point is approached with an exponent close to 2.2. In the 
critical gel they find algebraic decay of density fluctua- 
tions with an exponent value of « 0.35. All these findings 
agree qualitatively with the predictions of our model. In 
fact, our expression (||) for S t (q) has been suggested on 
phenomenological grounds as a starting point for a dis- 
cussion of the critical dynamics at the sol-gel transition 
. Although a general scaling of the diffusion constant 
of a individual cluster with its size was postulated, it has 
turned out that the so-called Rouse limit, which leads 
back to our representation (||), agrees well with experi- 
ment [ fl9f . In comparing our results with experiment one 
should keep in mind that we have used a distribution of 
crosslinks lacking any correlations, and have not taken 
into account excluded volume interactions. Presumably 
this amounts to a mean-field approximation of a short- 
range model. Vulcanization of a dense melt of linear 
macromolecules has been argued |2l|] to be mean-field- 
like, except within a narrow critical region of crosslink 
densities whose width decreases like the one third power 
of the ratio of the persistence length to the chain length. 
Since the argument includes and is in fact based on the 
distribution of cluster sizes, it should also hold for the 



dynamics. 
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